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We present constraints on testing general relativity (GR) at cosmological scales using recent data 
sets and assess the impact of galaxy intrinsic alignment in the CFHTLenS lensing data on those con¬ 
straints. We consider data from Planck temperature anisotropies, the galaxy power spectrum from 
WiggleZ survey, weak lensing tomography shear-shear cross correlations from the CFHTLenS sur¬ 
vey, Integrated Sachs Wolfe-galaxy cross correlations, and baryon acoustic oscillation data. We use 
three different parameterizations of modified gravity (MG), one that is binned in redshift and scale, 
a parameterization that evolves monotonically in scale but is binned in redshift, and a functional 
parameterization that evolves only in redshift. We present the results in terms of the MG parameters 
Q and E. We employ an intrinsic alignment model with an amplitude AcFHTLenS that is included 
in the parameter analysis. We find an improvement in the constraints on the MG parameters corre¬ 
sponding to 40 — 53% increase on the figure of merit compared to previous studies, and GR is found 
consistent with the data at the 95% confidence level. The bounds found on AcFHTLenS are sensitive 
to the MG parameterization used, and the correlations between AcFHTLenS and MG parameters 
are found to be weak to moderate. For all 3 MG parameterizations AcFHTLenS is found consistent 
with zero when the whole lensing sample is used, however, when using the optimized early-type 
galaxy sample a significantly non-zero AcFHTLenS is found for GR and the scale-independent MG 
parameterization. We find that the tensions observed in previous studies persist, and there is an 
indication that CMB data and lensing data prefer different values for MG parameters, particularly 
for the parameter E. The analysis of the confidence contours and probability distributions suggest 
that the bimodality found follows that of the known tension in the as parameter. 

PACS numbers: 95.36.+x,98.80.Es,98.62.Sb 


I. INTRODUCTION 

With the ongoing and future high precision surveys and missions such as the Dark Energy Survey (DES), Large- 
Synoptic Survey Telescope (LSST), Euclid and WFIRST, the question of testing general relativity at cosmological 
scales continues to drive a lot of interest in the quest to understand the nature of gravity and the dark energy associated 
with the observed cosmic acceleration. While probing gravity theories at cosmological scales is a legitimate endeavor 
by itself, it also tackles the question of whether the observed cosmic acceleration is due to an extension or modification 
to Einstein’s equations at cosmological scales, or is the result of a repulsive dark energy component permeating the 
universe. There exists a large body of literature on these questions. We provide a partial list references and reviews 
on the topic here HH3H, and we direct the reader to the citations therein. 

The precise results of cosmic microwave temperature anisotropies from Planck J3HHSD] , the recent data from the 
WiggleZ Dark Energy Survey[511 ;.52j, and the CFHTLenS weak lensing data [55H57] have all made it appealing to 
derive and analyze the latest constraints on general relativity or modified gravity parameters. These results will also 
help to guide us in advancing the analysis frameworks and to point us toward issues that need to be dealt with in a 
timely manner. 

In this paper, we contrast results from different ways to parameterize modifications to general relativity (GR). 
We generically call these modified gravity (MG) parameters, and we test several different types of parameteriza¬ 
tions including a binned parameterization that evolves in redshift and scale, a hybrid parameterization that evolves 


* Electronic address: jason.dossett@brera.inaf.it 

'Electronic address: 

tElectronic address: 

§ Electronic address: 


mishak@utdallas.edu 

d.parkinson@uq.edu.au 

tamarad@physics.uq.edu.au 



2 


monotonically in scale but is binned in redshift, and finally a functional parameterization that evolves only in redshift. 

In various analyses testing gravity at cosmological scales, weak lensing (cosmic shear) plays a prominent role in 
constraining the growth of large scale structure and thus the MG parameters. Weak lensing is also at the center 
of a number of future experiments aimed at testing gravity at cosmological scales (e.g. LSST, Euclid, WFIRST). 
However, an active area of work in the lensing community is also focused on understanding and controlling the 
systematic effects affecting this probe. At the forefront of these systematics are the intrinsic alignments of galaxies 
that generate correlations which contaminate the pure cosmic shear signal, see for example the reviews [581 [5!? ! 
and references therein. Briefly, there exist two types of galaxy intrinsic alignments. The first one is between close 
galaxies aligned with each other due the gravitational field present during their formation. These are referred to as 
intrinsic ellipticity - intrinsic ellipticity type or simply the II-type for the 2-point correlations and III for the 3-point 
correlations. The second type of intrinsic alignments are due to the fact that a massive structure aligns galaxies close 
to it and also produces lensing of background galaxies, resulting into an anti-correlation between cosmic shear and 
intrinsic ellipticities. This is known as the gravitational shear - intrinsic ellipticity type, or the Gl-type for the 2-point 
and can be generalized based on the same idea to the 3-point correlations giving the GGI and GII types. While 
the II and III intrinsic alignment can be suppressed by binning and cross-correlation techniques to assure that the 
galaxies are far enough and the gravitational tidal effect is small, these techniques cannot eliminate the GI, GGI and 
GII alignments since these are present between distant galaxies. Some theoretical analyses of the effect of intrinsic 
alignments and other lensing systematics on MG parameters can be found in [60[ Ifilj . We include in our analysis 
the possible effect of the 2-point II and GI intrinsic alignments present in the CFHTLenS data and analyze their 
correlations with the MG parameters in particular. 

There has also been increasing discussion about possible tensions between various cosmological data sets, partic¬ 
ularly in probing the amplitude of matter fluctuations at the CMB level versus probes at lower redshifts (e.g. weak 
lensing, galaxy clustering), see for example [B2] and references therein. The CMB tends to prefer models with higher 
values of the extrapolated cosmological parameter erg (the clustering amplitude on scales of 8h -1 Mpc) than that 
which is obtained from low redshift probes of the growth of structure. While some works have explored resolving 
these tensions with various changes to the neutrino sector (for example, |62H65j ). it is interesting to explore here if 
these tensions are present or reflected on the modified gravity parameters. 

The paper is organized as follows: in section II, we describe the methodology and parameterizations used. We 
describe the data sets used in the analysis in section III. The results and discussion are in section IV, while we 
conclude in section V. 


II. METHODOLOGY 

In order to constrain deviations from general relativity we update the publicly available package ISiTGR ESlHSj jijlT 
for use with the December 2013 version of CosmoMC m which is compatible with the likelihood codes for CMB power 
spectrum data from the Planck satellite [48]. Here we will briefly overview the modified growth formalism used in 
ISiTGR. A much more detailed account of the this formalism though is available in [251 ITS] . 


A. The Modified Growth Formalism of ISiTGR 

ISiTGR uses modified versions of the first order perturbed Einstein’s equations from the perturbed Friedmann- 
Lemaitre-Robertson-Walker (FLRW) metric. In a flat universe this metric is written in the conformal Newtonian 
gauge as: 


ds 2 = a(r) 2 [—(1 + 24')dr 2 + (1 — 2&)dx l dxi], (1) 

where r is conformal time, a(r) is the scale factor normalized to one today, the x^s are the comoving coordinates, 
and 4/ and $ are the potentials describing the scalar modes of the metric perturbations. 

Modified gravity parameters have been introduced in various interrelated notations that mainly parameterize a 
possible difference between the two potentials in the metric (for example, the gravitational slip parameter of 1331 ). 
and a second parameter that characterizes how the spacetime curvature side (or gravitational potentials side) is 
coupled to the source terms via the perturbed Einstein equations (sometimes this parameter can be related to a 
an effective gravitational constant). Due to degeneracies between some of these parameters, some other combined 
parameters have been proposed as exemplified at the end of this sub-section. We refer the reader to references [SH¬ 
UT]. a summary of the relationships between the various parameterizations can be found in [33] ■ It is also worth 
mentioning that papers have used these parameterizations in a functional or binned form. We use here the notation 
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Redshift bins 

Scale bins 

0.0 < z < 1 

1 < 2 < 2 

0.0 < k < 0.01 

Q 1, Ei 

Q 3 , E3 

0.01 < k < 00 

Q 2 , s 2 

Q4, S 4 


TABLE I: The subscript numbering in the binned parameterizations. 


introduced by m and used in ISiTGR 0B] as: 

/c 2 $ = —AnGa 2 '^^ piAi Q(k,a) (2) 

i 

k 2 (ty — R(k,a) $) = —12nGa 2 '^^p i (l + w i )aiQ(k,a). (3) 

i 

where, with i denoting a particular matter species, pi is the density, Aj is the rest-franre overdensity, and er, is the 
shear stress. Q(k,a) and R(k,a) are the time and scale dependent MG parameters, which both take a value of 1 in 
general relativity. Q quantifies a modification to what is often referred to as the Poisson equation (though as noted 
in [28] this equation is not truly the Poisson equation as it relates the overdensity to the space-like potential, $ which 
only affects relativistic particles). R then represents an inherent inequality between the two potentials that may be 
caused by a modified gravity model, and is known as the gravitational slip [551 . 

In order to not only avoid a strong parameter degeneracy between the parameters Q and R, but to also have a 
parameter that is directly probed by observations, ISiTGR does not directly use Eq. ^ in its code, but rather uses 
a combination of Eqs. © and © m- 

fc 2 (>F + $) = - 87 rGa 2 ^/ 9 ) A i E(fc, a) — 127rGa 2 pi(l + wi)<Ji Q(k, a), (4) 

i i 

where the parameter E = Q( 1 + R)/ 2. This parameter is directly probed by observations such as weak gravitational 
lensing. Like the parameters Q and R , E takes a value of 1 in general relativity. Note that in our previous works E 
was referred to as V, but in an effort to have a consistent set of parameters in the literature going forward we are 
now using the much more common, E. 


B. Evolution of the MG parameters 


In general there are two things that can be done with this modified growth formalism. First, one can give the 
MG parameters a generic form in order to look for possible deviations from general relativity. Alternatively, one 
can assume the MG parameters take a specific functional form in order to mimic the effects of a particular modified 
gravity model in order to test that particular modified gravity model under the assumption that the expansion history 
can be described by a wCDM model (a model where the dark energy has an equation of state, w). Since our goal in 
this work is to look for deviations from GR, we will of course be taking the first approach. 

We use three different MG parameterizations with different time and scale dependencies. We have described these 
parameterizations in our previous works [251 l46j and will briefly overview them again here. 

• PI: Firstly, we use a traditional binning parameterization in which the modified gravity parameters are binned 
in both redshift, z, and wavenumber (scale), k. A total of four bins are created by using two redshift bins and 
two scale bins. The scale bins are k < 0.01 and k > 0.01, while the redshift bins are 0 < z < 1 and 1 < z < 2. 
For redshifts z > 2 the MG parameters take their GR value of 1 at all scales. For continuity and numerical 
stability this parameterization is cast functionally as: 


with 


X{k,a) = ^(1 + X Zl {k)) + ^(X 2 (fc) - X Zl (k)) tanh^^ + ^(l - X Z2 {k)) tanh 


X Zl (k) 


X Z 2 (k) 


- (A 2 + Xl) + - (X 2 — Xi) tanh 

- (X 4 + X 3 ) + — (X 4 — X 3 ) tanh 


k- 0.01 
0.001 ’ 
fc-0.01 


(5) 


( 6 ) 


0.001 ’ 
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where X takes the values Q or E so in this parameterization a total of eight MG parameters are varied, E, and 

Qi, i = 1,2,3,4. 

• P2: Secondly, we use a hybrid parameterization that evolves with monotonic function in scale, but is binned in 
redshift identical to the P2, Eq. |5]). The scale dependence of this parameterization takes the form: 


X zi {k) = X-,e o.oi+x 2 (l-e rax), (7) 

X Z2 (k) = A' 3 e _ rar + X 4 (l - e _ rai), 

again giving a total of eight MG parameters, E, : and Qi, i = 1,2, 3,4. 

• P3: The third and final parameterization we use is what we have previously called the functional form param¬ 
eterization and was first introduced by m- This parameterization is scale independent and the parameters 
evolve only in time as: 


X(a) = {X 0 - 1 ) a s + 1 , ( 8 ) 

where X represents one of the MG parameters and Xq is the value of that parameter today. Our use of 
this parameterization though is slightly different from the two we discussed above in that we evolve the MG 
parameters Q and R with this functional form rather than Q and E. The evolution of E in this case is defined 
as described above using E(a) = Q(a)(l + R(a))/2. 


III. DATA SETS 

As a probe of the expansion history of the universe we use baryon acoustic oscillation (BAO) measurements from the 
6 dF Galaxy Survey measurement at 2 = 0.106 [58], the reanalyzed SDSS DR7 EM at effective redshift z e g = 0.35, 
and the BOSS DR9 [ 71 ] surveys at z e g = 0.2 and z e g = 0.35. These data sets help to break parameter degeneracies 
relating to late-time observables such as S 2 m . 

In order to access the wealth of information contained in the CMB, we use the measurements of CMB temperature 
anisotropy [45] from the first data release of the Planck surveyor. In the high-£ regime, the distribution of CMB angular 
power spectrum, Ci, can be well approximated by a Gaussian statistics, while the low-I? part of the CV distribution 
is non-Gaussian. For these reasons the Planck team divides the likelihood into low-£ (l < 50) and high-£ (£ > 50) 
parts and adopts different methodologies to build the likelihood in each region. The low -£ part of the likelihood 
employs a physically motivated Bayesian component separation technique to separate the cosmological CMB signal 
from diffuse Galactic foregrounds. All of the Planck frequency channels, from 30 to 353 GHz, are used for this part 
of the likelihood. The high- 1 7 part of the likelihood, however, utilizes a correlated Gaussian likelihood approximation 
based on a fine-grained set of angular cross-spectra derived from multiple detector power-spectrunr combinations 
between the 100, 143, and 217 GHz frequency channels. For this part of the likelihood, foregrounds are accounted for 
by marginalizing over power-spectrunr foreground templates. The first data release from Planck did not include any 
analysis of polarization data from the satellite’s observations, and for this reason the Planck likelihood code uses the 
low-1? WMAP polarization likelihood (WP) [43H7M74] . which is useful to break the well-known parameter degeneracy 
between the reionization optical depth r and the scalar index n s . Finally, unresolved foregrounds are marginalized 
over, assuming wide priors on the relevant nuisance parameters as described in [50] . 

As a first probe of large scale structure in the universe, we use measurements of the galaxy power spectrum as made 
by the WiggleZ Dark Energy Survey. The WiggleZ galaxy power spectrum is measured from spectroscopic redshifts 
of 170,352 blue emission line galaxies over a volume of 1 Gpc 3 EH Eg. The covariance matrices for the likelihood, 
as given in i 52] , are computed using the method described by ESI ■ To minimize the possible effects of non-linearities, 
we restrict ourselves to scales less that /c max = 0.2/i/Mpc and use the best model proposed for non-linear corrections 
to the matter power spectrum for this data set as calibrated against simulations (model G in [52]). Finally, we also 
marginalize over a linear galaxy bias for each of the four redshift bins, as in [52| . 

Our next probe of the structure of the universe is Integrated Sachs Wolfe (ISW)-galaxy cross correlations as 
presented in m- This data set cross correlates the CMB temperature anisotropies from the ISW effect with the 
galaxy distributions measured from the 2MASS and SDSS luminous red galaxy (LRG) surveys. It is a useful probe 
in that the galaxy distribution is sensitive to changes to the actual Poisson equation which probes changes in the 
parameter combination 2E — Q, while the ISW effect directly relates to the time derivative of the sum of the metric 
potentials and thus probes E and E. This likelihood and its usefulness in testing for deviations from general relativity 
has been discussed extensively in previous works, for example, [251 l4TT - [4il [46] . 
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As a final probe of the matter distribution of the universe we use the weak tensing tomography shear-shear cross 
correlation data. This data is not only useful in that it helps to constrain both the expansion history and the growth 
history of structure in the universe, but, as discussed above, it also gives us a way to directly probe the value of 
the MG parameter E. In this work we use tomographic shear cross-correlation data from the Canada France Hawaii 
Telescope Lensing Survey (CFHTLenS) [771 . The CFHTLenS survey analysis combined weak lensing data processing 
with THELI |53| . shear measurement with lensfit [55], and photometric redshift measurement with PSF-matched 
photometry [53i ■ Analysis of systematic and photometric redshift errors for the weak lensing data was discussed in 
[5B1EZ1- The data set of m consists of 21 sets of cosmic shear correlation functions associated with six redshift bins, 
each spanning the angular range of 1.5 < 9 < 35 arcrnin. The calculation of the theoretical shear cross correlation 
functions for this data set differ somewhat from those used in the weak lensing data set previously included in ISiTGR, 
m, so we will briefly overview the relevant equations for the calculation here. A more detailed description, though, 
can be found in m- 

As usual, the shear cross correlation functions £+ _(A)gg between bins k,l are given by 

1 r°° 

e + Uo) GG = — j dit j QA {ie)p*\t), (9) 

where J n is the n th -order Bessel function of the first kind, £ is the modulus of the two-dimensional wave vector, and 
P kl is the convergence cross-power spectra between bins k and l is given by m 

P *W = J d X9k(x)9l(x)P<l>,<t>{j^ ww.x), (1°) 

with comoving radial distance, x, comoving distance to the particle horizon, Xh, and comoving angular diameter 
distance, fr<(x)- Above, we have absorbed the usual extra terms into the power spectrum of the sum of the metric 
potentials, P^, where 4> = including the MG parameters. The weighted geometric lens-efficiency factor for the 

k th bin, 5 fc(x), is given by 


t \ - 1 f a > ( ~ X) 

9k{x)=-rX d XPk(X) f , a » (11) 

<HX) J x Mx ) 

corresponding to the normalized galaxy redshift distributions pk ■ 

In order to mediate the effects of possible intrinsic alignment contamination to the weak lensing signal, we follow 
the technique used in m which parameterizes the contribution of the intrinsic alignments to the shear correlation 
measurements using a nonlinear intrinsic alignment model introduced by [80] , which is based on the linear tidal field 
alignment model of [ST] , which in turn is based on earlier work of [|3]. 

In this technique the measured correlation functions are considered to be a sum of contributions from intrinsic 
alignments, which we will denote £+_(0 ) GI and £+ _(0 ) n , and the true shear correlation function, £+ _(0) GG above, 

£',-(*) = £+,—Wn + £+,-(0 )gi + £+,-Wgo- (12) 

The intrinsic alignment contributions, . _( 6 ») G i and ^_( 6 ») n , are calculated in the same way as £+_(0) GG , Eq. @ 
except that the convergence cross power spectrum, P K , is replaced by the projected GI and II power spectrum, P G1 
and P u respectively. Explicitly these power spectra are given by: 


(13) 


P kl it) = [ dx Pk ^x)P l (x) p2 p ( _ [ _ \ 

n l j Jo x ifM? 1 So ’ So ^fK(xy Xp 


\Mx)T F ' 

where do is the matter overdensity today and Pi is a cosmology dependent factor given by: 

Pi = — AcFHTLenS Cl Pcrit Hm ■ 


(14) 


(15) 


Above, p C rit is the critical density of the universe today, C\ is a constant with a value 5 x 10 14 h 2 M P) 1 Mpc 3 , and 
AcFHTLenS is a nuisance parameter that we will marginalize over in our likelihood analysis. 
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FIG. 1: 68% and 95% 2-D confidence contours for the parameters Qi and E, from parameterization PI for redshift and scale 
dependence of the MG parameters. All of the constraints for this evolution method are fully consistent with GR at the 68% 
level. 


IV. RESULTS AND ANALYSIS 

For all results we fit for the MG parameters and relevant nuisance parameters for the various data sets, as well as 
the six core cosmological parameters: and f l c h 2 , the baryon and cold dark matter physical density parameters, 

respectively; 9 , the ratio of the sound horizon to the angular diameter distance of the surface of last scattering; r, the 
reionization optical depth; n s , the spectral index; and lnl0 10 A s , the amplitude of the primordial power spectrum. 
As this analysis aims to test general relativity, primarily the ACDM model, we assume a ACDM expansion history 
throughout our analysis. This assumption should have little impact on our results. As shown in [S3] the MG parameters 
are robust to deviations from this expansion history, even those which include dark energy with perturbations. 


95% confidence limits on MG parameters 
evolved using form PI 

Qi 

[0.49,2.56] 

Si 

[0.97,1.14] 

q 2 

[0.05,3.08] 

£2 

[0.84,1.22] 

q 3 

[0.30,1.78] 

£3 

[0.97,1.06] 

Q4 

[0.28,2.88] 

£4 

[0.90,1.12] 


TABLE II: We list the 95% confidence limits for the MG parameters from using form PI to define their time and scale 
dependence. In this traditional binning approach we find that all of the MG parameters are fully consistent with their GR 
values of 1 at the 95% level. 


A. PI Traditional Binning 

We first, presents our results when using form PI for the evolution of the MG parameters. The one dimensional 
marginalized constraints for this parameterization are presented in Table [IT] We also show the 2 -D marginalized 68 % 
and 95% confidence contours for the parameters in Fig. [T} Using this evolution method, all of the MG parameters 
are consistent with their GR value of 1 at the 95% level and no noticeable tensions are evident. This is in contrast to 
the other evolution methods we use where some noticeable tensions appear, making those results more revealing as 
we will discuss more in the sections below. 


B. P2 Hybrid Evolution 

The results when using the hybrid evolution method for the MG parameters, P2, are more interesting. This can 
be seen quickly in Table |III| where we give the marginalized constraints on the MG parameters for this evolution 
method. Of particular interest are the constraints on the parameter Ej where the GR value of 1 lies outside of 
the 95% confidence interval. This is very interesting as it could signal a possible deviation from general relativity. 
Looking at the 2-D confidence contours for the MG parameters in Fig. [2] offers a bit of a reprieve though as the 
GR point still lies within the 95% confidence level in the 2-D parameter space. While it is well known that when 
using a large number of parameters, one of them may seem to indicate new physics when there is none (the so-called 
“look-elsewhere effect”), the tension exhibited in the constraints on this parameter, £i, which again is only one of 
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FIG. 2: 68% and 95% 2-D confidence contours for the parameters Qi and E. t from parameterization P2 for redshift and scale 
dependence of the MG parameters. As you can see in the first bin, there a tension with the GR value of 1. However, contrary 
to the marginalized 1-D constraints given in Table |Ih| the GR point is still within the 95% confidence region. 


eight MG parameters, should nevertheless be explored. We delay the majority of this discussion to §IVD[ where 
after presenting the results from evolution method P3 it becomes more clear that this constraint is driven by a long 
standing tension between the weak lensing and CMB data sets. 


95% confidence limits on MG parameters 
evolved using form P2 

Qi 

[0.38,3.43] 

Si 

[1.03,1.37] 

q 2 

[0.00,2.86] 

£2 

[0.75,1.07] 

q 3 

[0.28,2.46] 

£3 

[0.93,1.14] 

Q 4 

[0.05,1.99] 

£4 

[0.86,1.14] 


TABLE III: We list the 95% confidence limits for the MG parameters from using form P2 to define their time and scale 
dependence. While most of the MG parameters for this hybrid evolution method are consistent with their GR values of 1 at 
the 95% level, we find a significant tension for the MG parameter E in the large scale (small k ) low redshift bin (Ei). 


Correlation table 

Binning parameterization (PI) 


Q 1 

Q 2 

Q 3 

Q 4 

Si 

s 2 

S 3 

s 4 

AcFHTLenS 

-0.021162 

-0.29209 

0.015916 

0.0056355 

-0.0014863 

0.083586 

0.015755 

0.066954 

08 

-0.012168 

-0.53048 

0.044293 

-0.43088 

0.045781 

-0.61952 

0.048845 

-0.29894 


-0.0012586 

-0.072645 

-0.051569 

0.11762 

-0.08057 

-0.085185 

-0.033916 

-0.17292 

Hybrid parameterization (P2) 

AcFHTLenS 

0.058535 

-0.29535 

-0.052588 

0.095984 

-0.14858 

0.20636 

-0.086038 

0.10421 

08 

0.2655 

-0.70809 

0.12172 

-0.33026 

0.32713 

-0.59009 

0.1362 

-0.20504 


0.027229 

-0.065934 

-0.028016 

0.0803 

0.01565 

-0.15645 

0.14932 

-0.26513 


TABLE IV: Correlations between AcFHTLenSi °8 5 Pm versus the MG parameters of PI and P2. 


C. P3 Scale Independent Evolution 

The results from the scale independent evolution method, P3, are also much more revealing than those from PI. 
They offer some insight as to the origin of the noticeable tensions exhibited by parameters from P3. The marginalized 
95% confidence limits for the parameters Qm £o and Ro are presented in Table [V] The constraints are in general 
consistent with those of [ 2 S] , although a direct comparison is not very simple due to different additional probes used in 
the two works beside the CFHTLenS data. Now, looking only at these constraints, one might conclude that nothing 
interesting is happening in this parameterization of the MG parameters as all of the parameters are completely 
consistent with one at the 95% level. 

The 2-D confidence contours for these parameters plotted in Fig. [3j however, tell a very different story. Looking at 
these plots it becomes apparent that there is a noticeable tension in the parameter space, with the parameters showing 


























































95% confidence limits on MG parameters 
evolved using form P3 

Qo 

[0.77,1.99] | E 0 | [0.79,1.16] | R 0 | [-0.23,1.18] 


TABLE V: We list the 95% confidence limits for the parameters Q o, Eo, and Rq from the scale independent method, P3, 
of defining the evolution of the MG parameters. These constraints show no apparent tensions with GR. This is in contrast 
to the noticeable tensions that become evident when looking at the 2-D confidence contours shown in Fig. [3] The reason for 
this discrepancy is can be explained by looking at the 1-D probability distributions in Fig. [4] where the non-Gaussianity of the 
parameter constraints is easily seen. 



FIG. 3: 68% and 95% 2-D confidence contours for the parameters Qo, Eo, and R 0 from the scale independent parameterization, 
P3, for the MG parameters. These constraints are consistent with GR a the 95% level, but a tension is evident. The tension is 
evident when viewing these plots is not easily seen using the 1-D constraints given in Table [Y] This is due to the non-Gaussianity 
of the probability distribution for these parameters as further seen in Fig. [4] 


some preference for non-GR values. The constraints are also non-Gaussian. The non-Gaussianity of the parameter 
constraints is illustrated even better in Fig. [4] where we plot the 1-D probability distributions for the MG parameters 
used in this evolution method. The distributions for the parameters Q o and Rq are both quite skewed, with a large 
tail in one end of each of the distributions, while distribution for E 0 is almost bimodal. As discussed briefly in : |l V 1 1| 
the cause of these tensions in the MG parameter space seems to be a tension between the CMB and weak lensing 
data sets as can be seen in figure [5j We discuss this in depth in the next section. 


D. Tensions Between the CMB and Weak Lensing Data Sets 

As we have discussed above, there are noticeable tensions with GR in the constraints on the MG parameters from 
evolution methods P2 and P3. One of the eight parameters from P2 even shows an apparent deviation from GR at 
the 95% level when looking at the marginalized confidence limits. We have argued briefly above that this tension with 
GR in the MG parameter constraints arises from a tension between the CMB (Planck) and weak-lensing (CFHTLenS) 
data sets. We present a more in depth argument for that point here. 

There has been a known tension between CMB and weak-lensing data sets for quite some time (see and 
references therein). During an analysis of the data using the ACDM model, this tension is usually evident from 
comparing the preferred erg values of each of the data sets. While the CMB data usually prefers a erg ~ 0.83 [33] , 
lensing data usually prefers a much lower erg ~ 0.7. Recently there has been a lot of work on resolving this tension 
between these data sets by allowing for a sterile neutrino species or heavier set of active neutrinos (62H65j . In principle 
this would resolve the tension, as the neutrinos would suppress in late time growth and therefore explain why late 
time measurements of erg from lensing do not match those from the CMB, which infers the value of erg from early 
universe observations rather than directly measuring it. 

The biggest indication that the tension in the MG parameter space is coming from the known tension between the 
two data sets is the bimodal distribution of Eo from P3. This of course indicates that two very different values of 
Eo are equally preferred by the overall combination of data sets. Since erg has been useful to illustrate the tension 
between these two data sets before, we explore the values preferred when using P3. We find very interesting results 
upon doing this. 

In Fig. [5] we show the 2-D confidence contours in the Eo, erg plane as well as the 1-D probability distribution of erg 
and Eo for different combinations of the data sets that highlight the tensions between CMB and weak lensing (WL). 
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Correlation table 

Functional parameterization (P3) 


Q 0 

E 0 

AcFHTLenS 

-0.023164 

0.10624 

a s 

-0.66775 

-0.75738 


-0.072171 

0.052317 


TABLE VI: Correlations between AcFHTLenS, 08 , versus the MG parameters of P3. 





FIG. 4: 1-D probability distributions for the parameters Q o, Eo, and Ro from the scale independent parameterization for 

the MG parameters, P3. These constraints are consistent with GR a the 95% level, but a tension is evident particularly in 
the parameters Qo and Ro- The probability distributions display a significant level of non-Gaussianity, especially Eo which 
is somewhat bimodal. This non-Gaussianity explains why the tension seen in these plots is not seen using the marginalized 
constraints given in Table |Y] 


The first thing that can be noticed is as, even more so than Eo is bimodal. In fact the two peaks in the distribution 
correspond roughly to the preferred a 8 values of the CMB and weak lensing data sets. This is because, compared 
to GR, the MG parameters allow the CMB to fit the Planck data better with lower values of a$ and at the same 
time allow for a better fit to the the weak lensing data with higher values of a$- Quite importantly, from the 2-D 
confidence contours we see the higher values of Eo occur in the region of the parameter space where the lower a$ 
values are preferred. This lends credence to the conclusion that the higher preferred Ei from P2 is also being caused 
by this tension between the CMB and weak lensing data sets. 

While the tensions with GR we have seen in this work are still weak, they are nonetheless more significant than we 
have seen in the past. The fact that it is primarily due to a tension between two cosmological observations probing 
different eras of cosmic evolution is also important. Other explanations for the tensions between these data sets such 
as changes to the neutrino sector have been proposed, though some arguments have been made that these do not 
completely explain the tensions seen [521 . It remains an open question whether these observed tensions are signaling 
some issues with the observations or underlying theory. 


E. Galaxy Intrinsic Alignments: Constraints and Correlations 

Our results are presented in Figs. ©0 Tab. |IV[ Tab. [Vljand Tab. |VII| When the underlying theory is fixed to GR, 
we find constraints on the the amplitude of the intrinsic alignment model, ^CFHTLenS that are consistent with zero for 
the CFHTLenS full galaxy sample, the blue or the red samples (Figs. §0. But we find a clear non-zero AcpHTLenS 
parameter when we use the optimized-red galaxy sample of m in agreement with their results. This is expected 
since the optimized sample uses red-foreground and blue-background galaxies thus maximizing the signal to noise of 
the GI measurement for the red galaxies, as found. For GR, we find for this sample Ac FH TLenS = 3-54 ± 0.98 for 
the 68% limits. We find in general similar results to those of GR when the scale independent MG parameterization, 
P3, is used. However, when we use the binned and scale dependent MG parameterizations, PI and P2, we find that 
the zero ^CFHTLenS parameter is on the boundary line of the 95% confidence contours. This is possibly due to the 
larger parameter space in these two cases leading to larger contours. From the correlations tables, there are overall 
only weak to moderate correlations between the MG parameters and the intrinsic alignment amplitude parameter. 
It is found that Q 2 and E 2 (i.e. smaller scale and lower redsliift bin) are the MG parameters most correlated with 
^CFHTLenS- This is the bin where most of the CFHTLenS data resides as can be seen in figure 1 of m- We find no 
bimodality in the parameter AcpHTLenS but different distributions depending on the MG parameterization with P3 
being the closest to that of GR. We also find only a moderate relative change in the figure of merit (i.e. 1.3 — 5.3%) 
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FIG. 5: To illustrate that the tensions in the MG parameters may be arising from tensions between the Planck and CFHTLenS 
data sets, we plot the 68% and 95% 2-D confidence contours between Eo and as as well as the 1-D probability distribution 
for as and Eo when using method P3 for evolving the MG parameters. We plot three different data set combinations: CMB 
which includes all the data except CFHTLens; WL which includes CFHTLenS but not the CMB data; and CMB + WL 
which is the combination of all the data. The probability distribution for as is significantly bimodal with peaks corresponding 
to the preferred values of <78 for the Planck and CFHTLenS data sets. From the plots of the 2-D confidence contours, one can 
see that the bimodality of Eo seen in Fig. [4] is directly tied to that of as with the CFHTLenS data causing a preference for 
lower values of Eo and CMB data from Planck having a preference for higher values of Eo. 


for the MG parameters when AcFHTLenS is fixed versus when it is varied showing a moderate effect but this is likely 
to change for future high precision survey. 


V. CONCLUSION 

In this work we have placed constraints on deviations from general relativity using the modified growth formalism 
and a combination of the latest cosmological data sets including: observations of the CMB anisotropy power spectrum 
from the Planck satellite; observations of the galaxy power spectrum from the WiggleZ Dark Energy Survey; weak 
lensing tomography shear-shear cross correlations from the CFHTLenS survey; ISW-galaxy cross correlations; and 
BAO observations from 6dF, SDSS DR7, and BOSS DR9. We have also included in the analysis the effect of galaxy 
intrinsic alignment as a systematic effect in the weak lensing data. We used a model with a single nuisance parameter 
AcFHTLens to account for intrinsic alignments. We expressed the results on modified gravity in terms of the parameters 
Q and E. 

The constraints obtained using these latest data sets give a 40 — 53% improvement on the figure of merit for the MG 
parameters compared to previous studies [463 • GR is found to be consistent with observations according to all MG 
parameters when the 2-D marginalized 95% confidence contours are considered. We derived bounds and correlations 


Correlation table 

Binning parameterization (PI) 


crs 


AcFHTLenS 

0.38292 

0.052433 

Hybrid parameterization (P2) 

AcFHTLenS 

0.32715 

0.046858 

Functional parameterization (P3) 

AcFHTLenS 

0.12834 

0.089536 


TABLE VII: Correlations between the intrinsic alignment amplitude parameter AcFHTLenS and the amplitude of matter fluc¬ 
tuations as as well as the matter density parameter Q m , for the three MG parameterizations. 

of amplitude of the intrinsic alignment, AcFHTLenS- We included results for CFHTLenS galaxies split into early-type 
(red) and late-type (blue), and the optimized early-type sample (red-foreground-blue-background to maximize GI 
intrinsic alignments) of [ 77 ]. The amplitude is found consistent with zero for the whole galaxy sample when the 
theory is fixed to GR as well as when modified gravity is allowed. But a significantly (95% CL) non-zero AcFHTLenS is 
found when the optimized early-type sample is used for GR (i.e. Ac FH TLenS = 3.54 ±0.98 for the 68% CL, consistent 
and slightly improved limits with m- We obtain similar detection results when our functional scale-independent 
MG parameterization P3 is used. For our binned scale-dependent MG parameterizations (PI and P2), the bounds 
are larger and the zero AcFHTLenS is on the border line of the 95% confidence contours, most likely due to a larger 
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FIG. 6: 68% and 95% 2-D confidence contours for the intrinsic alignment amplitude parameter TlcFHTLenS and Il m . FIRST 
ROW: the theory is fixed to GR and the constraints obtained are in good agreement with those of EZI though improved due 
to more precise recent data. To the left are the results for the intrinsic alignment optimized red galaxy sample of ]jT[- The 
clear detection of non-zero ^4cFHTLenS is also in agreement with |77| . SECOND and THIRD ROWS: Similar constraints are 
presented but for the scale-dependent parameterizations P2 and P3 that model any deviation from GR. The bounds are larger 
but a zero AcFHTLenS parameter is practically on the the 95% CL boundary line for the optimized red galaxy sample. FOURTH 
ROW: results for the scale independent MG parameterization. The constraints are very similar to the GR case with a robust 
non-zero .AcFHTLenS ■ This may hint to the effect of scale dependence in the MG parameterizations 
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0.7 0.8 

°8 


0.270 0.285 0.300 0.315 0.330 



FIG. 7: LEFT and MIDDLE: comparative 68% and 95% 2-D confidence contours for the intrinsic alignment amplitude param¬ 
eter TlcFHTLenS versus as, and AcFHTLenS versus f2 m , for the 3 MG parameterizations. RIGHT: 1-D probability distributions 
for the parameter AcFHTLenS- TOP ROW: constraints for the full CFHTLenS galaxy sample. BOTTOM ROW: constraints 
for the optimized red-galaxy sample. 


parameter space in the binned cases. We find overall weak to moderate correlations between AcFHTLenS and MG 
parameters, the largest ones being with our parameters Q 2 and £ 2 . Only a relatively moderate change in the figure of 
merit (i.e. 1.25 — 5.30%) for the MG parameters is found when we include or not intrinsic alignments in the analysis 
but the effect is likely to be more pronounced for future high precision survey. 

For our MG parameterization P3, we obtained a bimodal probability distribution for the MG parameter £0 indi¬ 
cating the presence of some tensions between the data sets. We analyzed some aspects of these tensions and find them 
related or similar to what is observed with the amplitude of matter fluctuation parameter as- Indeed, we find that the 
data seems to prefer a lower value of the modified gravity parameter £ 0 i n parts of the parameter space where higher 
values of og are preferred, while areas of the parameter space that prefer a lower value of as favor a higher value of 
So- This finding is consistent with tensions in the MG parameter space from previous studies and also with tensions 
in as reported in other recent works. Further work is needed in order to investigate the source of these tensions with 
the incoming higher precision data and surveys. 
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